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Zero-point fluctuations of surface plasmon modes near the interface between metal and nonlinear 
dielectric are shown to produce a thin layer of altered dielectric constant near the interface. This 
effect may be sufficiently large to produce multiple metastable states of the surface plasmon vacuum. 



Optical bistability in various nonlinear optical systems 
is a topic of great current interest because of its poten- 
tial implementations in optical computing^. The use of 
surface electromagnetic waves such as surface plasmon- 
polaritons^ is known to reduce considerably the optical 
power levels that are necessary to achieve optical bista- 
bility effects^. Such a power reduction is possible be- 
cause of considerable enhancement of the local optical 
field near the metal interfaces. Very recently a possible 
observation of optical bistability effect at single-photon 
levels due to localized plasmon excitation in nonlinear 
optical material-filled nanoholcs in gold films has been 
reported^. The experimental findings were confirmed to 
some degree by theoretical calculations in&t This devel- 
opment together with other photon-blockade effects^ ob- 
served in nonlinear optical cavities at single-photon levels 
may open new ways of designing optical quantum com- 
puters. It justifies continued interest in nonlinear optical 
effects at extremely low light levels. 

As we continue to reduce number of real quanta of en- 
ergy in a surface plasmon system, and still observe non- 
linear optical effects, we may ask ourselves if the prop- 
erties of the nonlinear dielectric near the metal interface 
are somehow altered due to zero-point fluctuations in the 
vacuum state of the surface plasmon system. In other 
words, can we observe nonlinear optical behavior even 
without a single plasmon quantum? It appears that the 
answer to this question should be yes. Zero-point fluc- 
tuations of the surface plasmon field E vac must produce 
some average E% ac (z), which is very weak at large dis- 
tances z from the metal surface, but must grow exponen- 
tially towards the surface (similar to the behavior of the 
E- field in the surface plasmon wave). Let us assume that 
the dielectric near the metal surface exhibits third-order 
nonlincarity: 

e d = e^+4n X ^E\ (1) 

where and x^ 3 ' are the linear dielectric constant 
and the third order nonlinear susceptibility, respectively 
(here we consider the case of a central symmetric nonlin- 
ear material with x^ — 0). Far from the metal surface 
e = e^. However, near the metal we should see an al- 
tered layer of the material with the dielectric constant 
shifted by 4ttx^ E% ac . But how large and noticeable is 
this effect? 



Our results indicate that this effect may be surprisingly 
large even in materials with modest x^ nonlinearities. It 
may be sufficiently large to produce multiple metastable 
states of the surface plasmon vacuum, so one can talk 
about optical bistability of the vacuum itself, regardless 
of the number of real photons or plasmons in the system. 
Even though surprising, this effect is not unique or ex- 
tremely unusual. Metastable vacuums do occur in many 
field theories such as quantum chromodynamics, quan- 
tum gravity, etc. The origin of our Universe is believed to 
be due to a spontaneous decay of some metastable quan- 
tum vacuum state. Thus, in addition to optical comput- 
ing applications, experimental study of metastable sur- 
face plasmon vacuums may provide us with a useful toy 
model which could facilitate our understanding of some- 
what similar field-theoretical models of other quantum 
vacuums. 

Let us start by considering the zero-point energy of 
the surface plasmon vacuum of a square a x a region of a 
thin metal film with thickness d « a (Fig. 1(a)). Let us 
assume that the half-spaces above and below the film are 
filled with the same third-order nonlinear dielectric and 
consider the dispersion law of a surface plasmon (SP), 
which propagates along the metal-dielectric interface (the 
SP field decays exponentially both inside the metal and 
the dielectric). In such a case the dispersion law can be 
written asS 

,2 _ ^ 2 e d e rn (u;) 

c 2 e d + e m (u)±2e d e-^ [Z> 

where e m (w) is the frequency-dependent dielectric con- 
stant of the metal. If we assume that according to the 
Drude model e m = 1 — uj^/uj 2 (lu p is the plasma frequency 
of the lossless metal), and d is large, the dispersion law 
can be simplified as 

k2 = ^j^M (3) 

c z t d + e m (w) 

where e m (w) is real. This dispersion law is shown 
in Fig. 1(b) for the cases of metal- vacuum and metal- 
dielectric interfaces. It starts as a "light line" in the 
respective dielectric at low frequencies, and approaches 
asymptotically uo S p — ^/(l + td) 1 ^ 2 at large wave vec- 
tors. The latter frequency corresponds to the so-called 
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surface plasmon resonance. The dispersion law (2) also 
looks very similar to Fig. 1(b) in the general case of a 
lossy metal film if the metal film thickness d is small, 
so that the imaginary part of the term ±2ede~ kd com- 
pensates the imaginary part of e m . In such cases the 
plasmon wavevector k also diverges^ at some frequency, 
which is close to ui sp = U) p /(1 + e^) 1 / 2 . Since at every 
wavevector the SP dispersion law is located to the right 
of the " light line" , the SPs of the plane metal-dielectric 
interface are decoupled from the free-space photons due 
to the momentum conservation law. Nevertheless, very- 
short wavelength (in the 10-50 nm range) plasmons with 
the frequencies near the frequency of surface plasmon res- 
onance were observed in a number of near-field optical 
experiments^*^. 

The main contribution to the zero-point energy of the 
surface plasmon field comes from the frequency region 
near the surface plasmon resonance u) sp in which w does 
not depend on k, since the plasmon density of states is 
very large in this region. As a result, the total zero-point 
energy of the plasmon vacuum can be approximated as 



U vac « -T, k hld sp , (4) 

where summation has to be done over all the plas- 
mon modes k of the square region of the thin metal film 
under consideration. The surface plasmon eigenmodes 
of this square region are defined by the two-component 
wave vector (k x ,k y ) = ir/a x (n x ,n y ), where n x and n y 
are integer. Thus, if a of the order of 100-200 nm is 
selected, the quasistatic approximation used in equation 
(4) is well-justified. The summation over all possible sur- 
face plasmon wave vectors has to be cut off at some k max , 
which according toi£ is defined by the Landau damping 
and equals to the electron Fermi momentum kF so that 
I k max |~ &p = (67r 2 n e ) 1 / 3 (where n e is the electron den- 
sity). However, if we are less optimistic about the total 
zero-point energy, we may assume that k max = 27r/A m i„ 
corresponds to the shortest surface plasmon wavelength 
Amm observed in the experiment. 

Thus, the total number of plasmon modes on both the 
top and the bottom interfaces of the metal film is ap- 
proximately 4k max a 2 / ir 2 , and 
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U vac = —^h' jJ spk, m ax a (5) 

On the other hand, zero-point energy of the plasmon 
field in this region can be expressed as 



Uvac = — a 2 < l/k sp >, (6) 

where < E 2 p > is the average of the square of the 
electric field of the plasmon zero-point fluctuations, and 
< l/ksp >= 2/k max is the average decay length of the 



plasmon field away from the metal surface, which is ob- 
tained by averaging 1/fc over all the surface plasmon 
modes of the square film region. As a result, the elec- 
tric field of plasmon zero-point fluctuations can be found 

as 



< E % >= -^s P k 3 max (7) 

IT 

The shift of the dielectric constant of the nonlinear 
optical material near the metal surface equals to 

Ae d = 4ttx (3) < E 2 sp >= 16 X (3) ^ sp kf nax = Atu sp (8) 

Assuming rather cautious X m in — 30nm, we obtain 
Aed ~ I0 9 x^ esu -1 . Taking into account that such non- 
linear dielectrics as 3BCMU and 4BCMU polydiacety- 
lene materials have x^ as large as 10~ 10 esui 3 -, very 
large zero- field induced shifts Aed ~ 0.1 may be expected 
near the interface between metal and such nonlinear di- 
electrics (note that these materials were used in the ear- 
lier photon-blockade experiments^) . If less cautious up- 
per limit | kmax |~ kF is assumed, the resulting shift 
would be much larger (Aed ~ 10 14 x esu -1 ), which 
would make this effect noticeable even in regular nonlin- 
ear optical materials. 

The predicted nonlinear dielectric constant shifts are 
so large that the effects of bistability may be expected 
for the surface plasmon vacuum state. The origin of this 
effect would be roughly the same as the origin of bistabil- 
ity effects in regular surface plasmon propagation^. Using 
the constant A from equation (8) which is proportional 
to x' 3 ' , the frequency of surface plasmon resonance may 
be written as 



Wsp ~ (l + e „)Va - (i + e (i)+AV> 1/2 ( ' 

Let us initially assume that an d A are constants 
independent of the optical frequency (this is not the case 
for 4BCMU polydiacetylcnc in which rapidly oscil- 
lates in the visible range from 450 nm to 700 nm 13 ). Un- 
der such an assumption the frequency of surface plasmon 
resonance is given by the following cubic equation: 

Auj% + (e« + 1 K 2 p - lu 2 p = (10) 

The graphic solution of this cubic equation is presented 
in Fig. 2 for the cases of positive and negative x'- 3 ''. Posi- 
tive x' 3 ' increases the dielectric constant in the immedi- 
ate vicinity of the metal surface and shifts the frequency 
of the surface plasmon resonance downward because of 
the zero-point fluctuations. Negative x' 3 ^ has an op- 
posite effect on the surface plasmon resonance: the di- 
electric constant near the metal surface decreases, which 
leads to increase of the frequency of the surface plasmon 



3 



resonance. In addition, another meaningful root of equa- 
tion (10) may appear at even higher frequencies, if this 
root is below the plasma frequency of metal u> p . Analyti- 
cal solution of the cubic equation (10) indicates that this 
interesting behavior occurs near A 2 = 4(e' 1 ^ + l) 3 /27Wp 
and the two roots are located around 



^2,3 



( e (l) +1)1/2 



(11) 



Thus, two states are possible for the surface plasmon 
vacuum when \^ — const < 0. The upper vacuum 
state is metastable and once created, it should eventu- 
ally spontaneously decay into the lower-energy vacuum 
state. Since x °f such materials as 4BCMU polidy- 
acetylene is not actually constant as a function of fre- 
quency, but rather quickly oscillates in magnitude and 
sign in the visible frequency ranged, other scenarios of 
bistable or even multi-stable plasmon vacuums are pos- 
sible. If ( w ) changes sign from positive to negative in 



the vicinity of u 



/(l + e* 1 )) 1 / 2 , the lower-energy 



vacuum state from the pair obtained for x^ < 0, which 
is situated above u>ip becomes metastable, and would 
eventually decay into a stable state located below wjp at 
X (3) > 0. 

The metastable states of surface plasmon vacuum may 
be created using the effect of bistability observed exper- 
imentally in the excitation of surface plasmons over an 
interface between metal and nonlinear dielectric 4 . Real 
physical surface plasmons excited in the upper energy 
state would create a necessary seed layer with an altered 
dielectric constant near the metal interface. After all the 
real plasmons decay or radiate (due to the losses in metal 
and surface roughness), the surface plasmon vacuum may 
stay in the metastable state for a while. Exactly how long 
this metastable vacuum state is going to last would be 
best addressed by a pump-probe type of optical exper- 
iment. Theoretical answer to this question would meet 
problems, which are somewhat similar to the problems 
in quantum gravity, neither of which has been solved so 
far. In order to understand this difficulty let us recall the 
well-known analogy between media electrodynamics and 
gravitation theory. 

The realization that solid-state toy models may help 
in an understanding of electromagnetic phenomena in 
curved space-time has led to considerable recent effort 
in developing toy models of electromagnetioi 4 * 1 ^ and 
sonioH^L 1 ^* 1 ^ black holes. In the case of media elec- 
trodynamics this may be possible because of an analogy 
between the propagation of light in matter and in curved 
space-times: it is well known that Maxwell equations 
in a general curved space-time background gij{x,t) are 
equivalent to the phenomenological Maxwell equations in 
the presence of a matter background with nontrivial di- 
electric and magnetic permeability tensors €ij(x,t) and 



;{x,t) = fiij(x,t) 



m{x,t) 



(12) 



In this analogy, some features of the black hole event 
horizon may be emulated by a surface of singular di- 
electric constant and magnetic permeability, so that the 
speed of light goes to zero, and light is "frozen" near such 
a surface. 

If we accept this language, the nonlinear optical cor- 
rections to 6ij(x,t) and Hij(x,t) of the form described 
by eq.(l) would describe an effective gravitational in- 
teraction between the quanta of electromagnetic field. 
This fact is very easy to illustrate in situations in which 
€ij(x,t) and Hij{x,t) depend only on one spatial coordi- 
nate z (such as the Rindler geometry). After straight- 
forward calculations one obtains the following expression 
for the Rqq component of the Riemann tensor (recall that 
the trace of the energy-momentum tensor T = for the 
electromagnetic field): 



i d 2 i_ W% 2 2 



87rfc, 



eff 



(13) 



which reproduces the Einstein equation with the ef- 
fective gravitational constant k e f f , which is proportional 
to x^ ■ Unfortunately, in normal three-dimensional op- 
tics situations and x^ of optical materials are small, 
and no non-trivial effects of effective gravity may be em- 
ulated. All the electromagnetic black hole models con- 
sider motion of light in some static geometry of Uj{x) 
and fj,ij(x). 

Very recently we have pointed out that this situa- 
tion changes in the case of two-dimensional optics of 
surface plasmon-polaritonsSi. It is clear from Fig. 1(b) 
that around the frequency of surface plasmon resonance 
both phase and group velocity of surface plasmons tend 
to zero, and the effective two-dimensional dielectric con- 
stant of the dielectric diverges as seen by the plas- 
mons which propagate along the metal-dielectric inter- 
face. In a similar fashion the effective two-dimensional 

(3) 

X2 °f the nonlinear dielectric diverges near the sur- 
face plasmon resonance: small changes of the three- 
dimensional dielectric constant td due to the Attx^E 2 
term in equation (1) are perceived as very large changes 
of the effective two-dimensional dielectric constant by 
surface plasmon-polaritons. Thus, effective gravity of 
surface plasmons becomes very strong. As has been de- 
scribed in2£, gravitational collapse of surface plasmon 
field can be emulated in a self-focusing experiment per- 
formed near the frequency of the surface plasmon reso- 
nance. However, such self-focusing of surface plasmons 
emulates only the effects of classic general relativity 

In order to emulate quantum gravity one would need to 
make plasmons behave as the quanta of the effective grav- 
itational field. This is easy to achieve in a non-symmetric 
nonlinear optical medium with non-zero x 2 susceptibility. 
In such media 
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e d = e (1) +4ir X {2) E + 4ir X (3) E 2 1 (14) 

where the tensor indices of e, X 2 , X 3 > an d ^ are omit- 
ted for the sake of simplicity. Comparison of equations 
(12) and (14) indicates that oscillations of the effective 
gravitational metric are proportional to the oscillations 
of the E field of the surface plasmon waves. In effect, 
surface plasmons become effective gravitons of the sys- 
tem. Since the inversion symmetry is broken near every 
interface, surface x 2 terms are bound to appear at the 
interface of metal with any nonlinear optical medium, 
so that the equation (14) correctly describes nonlinear 
behavior of any such interface. 

To complete our quantum gravity analogy let us derive 
the effective Planck scale L e J^ = (Kk e f / /c 3 ) 1//2 of our 
system. Taking into account the value of the effective 
gravitational constant k e ff from equation (13), we find 

reff _ Meff 1/2 _ 877^X^2 

L pl ~ [ ~^3~> ~ ^(1)3 > [li> > 



Assuming x^ ~ 10 10 esu, and the plasmon wave- 
length A = 30 nm, which defines k — 2-ir/X, we obtain the 
effective Planck scale as L e JJ = 2 nm, while at A = 10 nm 

we get L e JJ w A. These simple arguments indicate that 
the effects of effective quantum gravity are very strong in 
our system. Since the lifetime of the metastable surface 
plasmon vacuum is determined by the quantum fluctu- 
ations of the effective metric of the nonlinear opti- 
cal material, theoretical resolution of this problem seems 
to be difficult at the moment. On the other hand, this 
lifetime may be measured experimentally in a toy "big 
bang" -like experiment. Such experiments would be very 
beneficial both from the point of view of designing opti- 
cal computers, and from the point of view of advancing 
theoretical models of quantum gravity. 
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FIG. 1: (a) Model geometry for the calculations of zero-point 
energy of surface plasmon vacuum in the case of a square ax a 
region of a thin metal film with thickness d « a. The half- 
spaces above and below the film are filled with the third-order 
nonlinear dielectric, (b) Surface plasmon dispersion law for 
the cases of metal- vacuum and metal-dielectric interfaces. At 
large plasmon wavevectors k the surface plasmon frequency 
does not depend on k. Nonlinear shifts of the surface plasmon 
resonance are shown for the cases of x^ > and X ^ < 0. 



FIG. 2: Graphic solution of equation (10) which defines the 
frequency of surface plasmon resonance for the cases of posi- 
tive and negative X ^ of the nonlinear optical material. The 
roots are shown by black dots. If x^ changes sign near lj 3P all 
roots become metastable states of surface plasmon vacuum. 
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